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In a recent papjsr an implicit equation for contacting viscoelastic 
spheres was derived Integrating this equation it can be shown that 
the coefficient of normal restitution e depends on the impact velocity 
g as 1 — e ~ 

The behavior of granular gases has been of large scientific interest in recent 
time. Goldhirsch and Zanetti □ and McNamara and Young 13 have shown that a 
homogeneous granular gas is unstable. After some time one observes dense regions 
(clusters) and voids. To evaluate the loss of mechanical energy due to collisions one 
introduces the coefficient of (normal) restitution 

9' = (1) 

describing the loss of relative normal velocity g' of a pair of colliding particles after 
the collision with respect-to the impact velocity g 

In the investigations □'cl the approximation of constant coefficient of restitution 
was assumed. Solving viscoelastic equations for spheres currently it was shown that 
the coefficient of xYoxsnijl restitution e is not a constant but a function of the impact 
velocity e{g) itselfBa. For the "compression" ^ = -Ri + i?2 — |n_ — ^2! of particles 
with radii i?i and R2 at positions ri and r2 one finds 

i + p(e" + lA^U)-0 (2) 



3 mcff (1 - v'^) 



mim2 ^ R1R2 

WefT = ■ -Keff = 



(3) 



TOi + TO2 i?l + i?2 
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Y is the Young modulus and m^fj and -R^ff are the effective radius and mass of the 
grains, respectively, yl is a material constant depending on the Young modulus, the 
viscous constants and the Poisson ratio of the material. The initial conditions for 
solving (||) are ^(0) = and ^(0) = g. The coefficient of restitution e of at time 
t = colliding spherical grains can be found from this equation relating the relative 
normal velocities g — f (0) at time of impact and at time tc, when the particles 
separate after the collision, i.e. tc is the collision time: 

e = -Htc)/iiO). (4) 

The (numerical) integration of equation (^) yield^i the coefficient of restitution as a 
function of the impact velocity (see fig. 1 in ref. El) which is in good agreement with 
experimental dataJ. Constaiit coefficient of restitution, however, does not agree 
with cjpjyfrimental experience U. Other theoretical work on this topic can be found 
e.g. in 



The duration of collision t^ for the undamped problem {A — 0) is given by □ 

©c 

— 2 1 I (5) 

with 8^ being a constant. Substituting ^ — p^^x{Q) we get for the rescaled ve- 
locity V = p^g, and with (^) one finds t — Ou^s. We rewrite @) using rescaled 
compression, relative velocity and time a;, v and O and the abbreviation a — 

X + av ^x\/x + v ^x^ —0 (6) 
with X = -^x. The initial conditions read now 

x{0) = (7) 

— (0) = ui — (0) = w, hence i;(0) = (8) 

Eq. (^ will be solved now by series expansion. All derivatives of third order and 
higher of x diverge at = 0, hence one cannot expand x in powers of 0. Since ^ 
can be written in the form ^(i) — gty(t), the Ansatz 

a;(e) = w^e (1 + 77(6)) with 77(0) = (9) 
seems to be a reasonable assumption. We get an equation for 77 : 

efj + 2ri + av^e^yjl + 77 + (av^VO + Q^'j (1 + t;)^ = . (10) 

In (|o|) occur terms 0°-^ and Q^'^, therefore we expand 77 in powers of \/Q 

00 

r/^J^a^e^ (11) 

fe=0 

The first coefficients oq = and ai = vanish because of the initial conditions. 
With Taylor expansion of ^/T+Tj and (1 + 7^)2 for small 77 we arrive at 
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■^avie^ - -^e^ + -^av^Q^ + -^o^viQ^ . . . (12) 

and therefore 



^ 15"" 35" ■ 70"- ' 15~ 



i ^ 4 5 4 4 7 1 2 ^r^A 

X = u5 6 - —av 62 - — ti5 92 -I a v'^Q 

15 35 15 

+la^;e5_ JLa3„le^ + — t;te6 + ... (i3) 

70 2475 175 ^ ' 

Rearranging the fuU series (|l3| ) one finds (s. fig. 0) 

a: = t)5xo(6) + avxi{Q) + a^u§2;2(9) . . . (14) 

v^xq is the sohition of the undamped (elastic) colhsion. We wiU need x (i0") 
where is the duration of the undamped colhsion. 
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FIG. 1.: The dynamics of the colhsion. The dashed Une shows the (strictly symmetric) solution of 
the undamped collision. For the case of the damped motion (full line) the maximum penetration 
depth is achieved earlier whereas the duration of the collision is longer (0c > ©2). The figure 
gives a nice impress why we employ the inverse collision instead of direct calculation: one cannot 
expand x beyond 0^. In the point 0^! the derivatives of the curve for elastic motion diverge. 



Later we will need the solution x™" of the inverse problem, i.e. of a collision with 
impact velocity v' and final velocity v. Hence the inverse collision is not a damped 
motion, but an accelerated one, where a has to be replaced by —a. Substituting 
V v' and a —a we find 

x™"(e') = (w')^xo(e') - av'xiie') + a''{v')ix2ie') . . . (15) 

Now we determine the collision time 8c and the final velocity -^x{Qc). (A more 
direct method to calculate Qc is to determine the solution of x{<d) — using Taylor 
expansion of x in the region close to 9°. It can be seen easily that this method fails 
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since, all derivatives of with n > 3 diverge for 8 = . Therefore 8c has to be 
calculated indirectly as explained in the text (s. fig. |l]).) 

The problem will be subdivided into two parts (s. fig. a) the motion of the 
particles x from 8 = to time 8m when x approaches its maximum and where x 
changes its sign, and b) from 8m to 8c. In case of undamped motion where a = 
we have 8m = 8^/2. In part b) we do not consider the collision itself but the 
inverse problem in the interval (8 = 0, 8m), with 8m being the time where a;*"" 
approaches its maximum. The continuity of both parts means 

x(8m) = x-(e:„) . (16) 




FIG. 2.: Sketch of the calculation. The first part (a) © G (0,0m) is calculated directly, for the 
other part (6) we define the inverse collision where the particles start with velocity v' and velocity 
approaches zero at = ©m. Both curves have to fit together smoothly. 

For finite damping a 7^ we write 8m = 8°/2 + (5 and 8;„ = (8°)' /2 + 6' and 
remind that = (8^)' . To get an expression for S we expand 

... /8"\ 6^ /80 



.M^o(-j+^^o^-j+-^.o^-j+...) (18) 
+va { xi { ^ ] + Sxi { ^ ] + —-^xi ^ 

(19) 



2/ V2/2d83V2 



and using xq (8c /2) = {v^ xq is the solution of the undamped problem) 

j ii (eg/2) 

xo (80/2 



' = ---^T:T^ + <^i-') ■ (20) 
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The expression (^0|) has to be inserted into the Taylor expansion of x (6c/2 
- + 5)^v*: (.0 (f ) + ) + (f ) + . . .) 



2 / V 2 / 2 io (eo/2) 



Hence 



a; (e™) = f = a:o (e°/2) + avxi (e°/2) + 

'e|\ 1 x[_(e2/2) \ 

2 ; 2 xo (90/2) j 



2 £ / 



Replacing again v ^ v' and a ^ —a yields 
„i i^i(9g/2) 

X- (e:„) = {v')^^Xo (e°/2) - (e°/2) + 



2 ; 2 io (02/2) y ' ■ ■ ■ 

As explained above both solutions (ESh and (E3) have to be equal. With 



^ "^^Uy 2x0 (90/2) 

we write 

t'^xo ( 1 +avxi ( ^ 1 +a^w5/3 = (w')5xo ( ^ 1 - av' xi ( 

+ a2(w')^/3 ■ 

We expand v' in a 



w' = w + avi + q:^V2 + 



and find 



with 6v = avi + a^V2 4 
coefficients yields finally 



Writing (1 + in powers of a and comparing 



f + —av^ 



f5 



y " ^0 (^) 

V {1- av^Ci + a^viC2 + 



H a v''- 
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and for coefficient of normal restitution one gets 



Xl 



m 



Xq 



(30) 



— = 1 - aw 5 Ci + a^w 5 C2 + . 

V 

3^.21 9^.942 



(31) 
(32) 



with g being the impact velocity (fig ^) . For the duration of collision we find with 
■"), (H) and m 



= e'^p-ig -t f 



1 

10 



Ciap^-g'i 



(33) 
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FIG. 3.: The coefficient of restitution over impact velocity due to eq. (|32|). As expected for small 
relative velocity the particles collide almost elastically. The result of numerical integration of (U) 
coincides with the curve. Both curves cannot be distinguished in the plot. 



Our final results eq. ( p^ ) shows that for viscoelastic colliding smooth bodies the 
coefficient of normal restitution is a decreasing function with rising impact velocity: 
1 - e ~ 

The authors are grateful to N. Brilliantov, S. Esipov, H. Herrmann, F. Spahn 
and W. Young for helpful discussions, and J.-M. Hertzsch for providing relevant 
literature. 
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